function solution=FE_solver_1D_Poisson(?)
%basis_type_trial==101: 1D linear

[P,T] = generate_PT(?,101);
if basis_type_trial==101
    
    Pb=?;
    Tb=?;
 
elseif basis_type_trial==102
    
    ?;
    
end



boundarynodes = generate_boundarynodes(?);



A = assemble_matrix_1D(?);


b = asemble_vector_1D(?);



[A,b] = treat_Dirichlet_boundary(?);


soultion=A\b;


error=compute_max_FE_nodes(?);